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I  Introduction 


A  hybrid  combination  of  the  geometrical  theory  of  diffraction  (GTD)  [1] 
-[4]  and  modal  methods  are  employed  within  the  framework  of  the  gen¬ 
eralized  scattering  matrix  technique  (GSMT)  [5]  to  efficiently  analyze  the 
electromagnetic  scattering  by  an  open-ended,  dielectric  filled  rectangular 
cavity  in  a  perfectly  conducting  ground  plane  where  the  medium  external 
to  the  cavity  is  assumed  to  be  free  space.  This  analysis  is  described  here  for 
the  two  dimensional  (2-D)  cavity  of  Figure  1  in  free  space  which  is  excited 
by  either  a  TE^  or  TM^  plane  wave;  however,  the  concepts  involved  can  be 
readily  extended  to  treat  three  dimensional  rectangular,  circular  or  other 
separable  cavity  configurations  by  augmenting  the  GTD  with  other  asymp- 
lotic  high  frequency  methods  as  shown  in  (6j.  An  analysis  of  the  latter 
three  dimensional  (3-D)  case  will  be  reported  separately.  Even  though  the 
width  of  the  cavity  is  shown  to  be  much  larger  tiian  its  depth  in  Figure  1, 
this  is  not  a  restriction  in  the  present  analysis  which  can  also  accommodate 
cavity  depths  that  are  larger  than  its  width. 

In  this  hybrid  analysis,  the  elements  of  the  generalized  scattering  ma¬ 
trices  associated  with  the  junctions  at  the  open  front  end  and  the  back 
termination  are  found  in  a  conceptually  simple  manner  via  liigh  frequency 
techniques  based  on  the  GTD  and  the  equivalent  ctirrent  method  (ECM) 
[7,8|.  It  is  noted  here  that  the  ECM  employs  the  GTD  to  calculate  Ihe 
strengths  of  the  equivalent  currents.  It  is  further  noted  that  the  high  fre¬ 
quency  techniques  can  be  used  because  each  waveguide  cavity  mode  can  be 
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decomposed  into  a  set  of  modal  rays  [9|.  The  interaction  of  these  waveguide 
modes  with  the  discontinuitues,  such  as  the  open  end,  of  the  cavity  can  then 
be  analyzed  via  GTD  and  ECM  by  using  the  diffraction  or  launching  of  the 
corresponding  modal  rays.  The  elements  of  the  scattering  matrices  furnish 
the  complex  amplitudes  of  the  waveguide  cavity  modes  which  are  reflected 
and  transmitted  by  the  junctions.  The  multiple  interactions  of  the  modal 
waves  between  the  junctions  are  fully  accounted  for  by  the  self-consistent 
multiple  scattering  approach  which  is  essentially  contained  in  the  GSMT. 
It  is  important  to  note  that  the  expressions  for  the  scattering  matrices  ob¬ 
tained  via  the  GTD  are  in  a  relatively  simple  form;  it  is  for  this  reason 
that  the  present  hybrid  approach  is  very  efficient.  In  contrast,  the  use  of 
conventional  mode  matching  procedures  to  analyze  this  problem  would  be 
far  more  cumbersome.  Also,  the  present  hybrid  method  of  solution  is  also 
far  more  efficient  than  a  moment  method  solution  of  this  problem.  On  the 
other  hand,  the  present  analysis  is  approximate  because  GTD  is  an  asymp¬ 
totic  method;  however,  it  is  also  quite  accurate  as  the  numerical  results 
based  on  this  approach  are  shown  to  agree  well  with  those  based  on  an  in¬ 
dependent  moment  method  solution  for  this  problem.  At  the  present  time, 
the  directions  of  incidence  and  scattering  are  restricted  to  about  10°  away 
from  grazing  with  respect  to  the  ground  plane;  further  study  to  remove  this 
restriction  will  be  pursued  in  the  future. 

The  format  of  this  paper  is  as  follows.  First  a  self-consistent  formulation 
of  the  GSMT  is  presented  in  section  II.  This  self-consistent  approach  [10] 
for  summing  the  multiple  wave  interactions  between  the  junctions  formed 
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by  the  open  front  end  and  the  back  termination,  respectively,  is  slightly 
different  from  the  Neumann  series  summation  approach  used  in  the  original 
development  of  the  GSMT  [5];  however,  it  leads  to  the  same  result  for 
the  multiple  interactions  as  in  the  original  GSMT.  The  various  junction 
scattering  matrices  are  also  defined  in  section  II.  In  section  III,  the  GTD 
and  ECM  are  employed  to  calculate  the  elements  of  the  junction  scattering 
matrices.  The  scattering  matrices  are  combined  in  section  IV  according 
to  the  expressions  for  the  multiple  wave  interaction  of  section  II  to  obtain 
numerical  results  for  the  fields  scattered  by  the  cavity  when  it  is  illuminated 
by  a  plane  wave.  Numerical  results  are  presented  for  the  backscattering  as 
well  as  for  the  bistatic  scattering  from  the  cavity  shown  in  Figure  1.  An 
time  convention  is  assumed  and  suppressed  throughout  this  analysis. 

II  Self-Consistent  Formulation  of  the  GSMT 

In  the  GSMT  [5],  the  basic  scattering  mechanisms  for  the  configuration  in 
Figure  1  can  be  isolated  and  identified  as  being  associated  with  the  junctions 
formed  at  the  open  front  end  and  the  back  termination,  respectively.  At 
the  present  time  the  GTD  edge  diffraction  coefficients  are  available  for  a 
conducting  wedge;  however,  the  edges  at  the  open  end  in  Figure  1  are  also 
in  contact  with  the  dielectric  which  fills  the  cavity.  Thus,  an  additional 
junction  is  artificially  introduced  in  the  ansdysis  so  that  the  GTD  for  a 
conducting  wedge  can  be  used;  this  artificial  junction  can  be  formed  by 
recessing  the  dielectric  by  a  small  distance  6  inside  the  cavity  as  shown  in 
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Figure  2.  An  analysis  of  the  problems  in  Figure  2  leads  to  a  solution  for 
the  scattered  field  of  the  problem  in  Figure  1  as  ^  >  0;  such  a  technique  of 

introducing  a  6  which  is  then  Edlowed  to  become  vanishingly  small  has  been 
employed  successfully  in  the  past  to  analyze  many  scattering  problems  via 
the  GSMT  [5,11]-  Three  junctions  are  now  identified  at  the  discontinuities 
as  shown  in  Figure  2;  these  junctions  are  labeled  I,  II  and  III,  respectively. 

The  scattering  matrices  associated  with  the  junctions  I,  II  and  III  de¬ 
scribe  the  isolated  scattering  events  at  these  junctions  as  shown  in  Figure 
3.  The  GSMT  extends  the  concept  of  the  conventional  scattering  matrix 
for  propagating  waveguide  modes  to  include  evanescent  modes,  and  it  also 
includes  the  exterior  scattering  matrix.  The  various  scattering  matrices 
associated  with  junctions  I  and  II  are  shown  below;  it  is  noted  that  in  this 
special  case,  the  effect  of  junction  III  is  combined  into  the  scattering  matrix 
characterizing  junction  II  and  is  hence  not  treated  separately. 

The  multiple  scattering  between  the  junctions  is  calculated  via  a  self- 
consistent  procedure  which  requires  a  knowledge  of  the  generalized  scatter¬ 
ing  matrices  [5n],  [5i2],  [‘S'2i],  [522],  and  [•S'rj  for  the  junctions  I  and  II. 

Let  E'  denote  the  electric  field  of  the  incident  electromagnetic  plane 
wave.  The  electric  field  contains  only  the  z  component  for  the  TM^  case 
and  it  contains  only  the  6  component  for  the  TE^  case.  Thus,  one  may 
express  E'  in  general  as; 


=  EU  +  Eii-ff) 


(1) 
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with 


Ei  = 

(2) 

E'e  = 

(3) 

where 

Je'  =  —ko{cosd'x  +  sin6'y),  (4) 

and 

p  =  xx  +  yy.  (5) 

It  is  noted  that  ko  refers  to  the  free  space  wavenumber  and  p  is  the  position 
vector  of  the  observation  point  P.  The  unit  vector  9'  is  shown  in  Figure  1 . 

The  matrix  [5ji]  relates  the  scattered  electric  field  E'°  from  only  the 
open  end  at  i  =  0  to  the  field  E'  which  is  incident  at  the  open  end.  One 
may  then  describe  the  scattering  matrix  [5ii]  as  follows: 


where 

=  El^z  +  E;°9  (7) 
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and  the  unit  vector  0  is  also  shown  in  Figure  1.  From  (6), 


ao  Cto 


(8) 


The  scattering  matrix  [5i2]  converts  the  waveguide  modal  fields  incident 
at  the  open  end  (i  =  0)  from  the  interior  region  to  the  field  radiated  by 
these  modes  from  the  open  end  as  shown  in  Figure  3(b).  The  modal  electric 
field  within  the  waveguide  portion  of  the  cavity  may  be  represented  as 


=  E  +  ^net)  (9) 

n  —  Q 

The  terms  and  are  the  amplitudes  of  the  modal  fields  in  the 
waveguide  section.  The  denotes  the  n*'*  modal  electric  field  in  the  TE, 
case  for  which  the  magnetic  field  had  only  a  i-component.  Likewise,  e^z 
denotes  the  modal  electric  field  which  is  entirely  in  the  z-direction  for 
the  TMj  case.  The  propagation  constants  /?„  for  the  mode  are  the  same 
for  the  TEj  and  TM^  cases.  Also  note  that  n  >  0  for  the  TM^  case.  The 
superscripts  +  and  —  in  (9)  refer  to  modes  propagating  in  the  +x  and  —x 
directions,  respectively.  These  waveguide  modes  are  defined  in  detail  in 
Appendix  A.  If  E^°  denotes  the  electric  field  at  P  exterior  to  the  cavity 
due  to  the  radiation  of  the  modes  that  impinge  on  the  open  end,  then  the 
scattering  matrix  [5i2]  relates  E'‘°  to  E^  as  follows: 
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(10) 


Ero 

■  [sn.l 

[s:i.)  ■ 

e;o 

.  [SJ3.1 

.  f^n]  . 

where 


[S-B-l 

f  1 

r  Cro  ] 

(11) 


Notice  that  each  of  the  submatrices  in  (11)  separately  contain  the  z  and 
6  components  of  the  fields  radiated  by  each  of  the  TM^  modes  and  the 
TEj  modes  respectively.  For  example,  each  element  of  gives  the 

i  component  of  the  electric  field  radiated  by  each  of  the  TM*  mode  of 
amplitude  Therefore,  if  the  total  numbers  of  the  TM^  and  TE^  modes 
included  in  the  calculation  are  N;,  and  N,,  respectively.  The  order  of  the 
submatrices  [Ssa„]^  !-®n]  will  be  IxN^, 

IxN/,,  IxN,,  IxN,,  Nf,xl,  and  N,xl,  respectively. 

Similarly,  the  scattering  matrix  [52i]  describes  the  coupling  of  the  inci¬ 
dent  plane  wave  field  into  the  waveguide  modes  as  shown  in  Figure  3(b). 
Therefore,  [52i]  relates  E~  to  E'  as  follows: 


'  [^n-]  ' 

[Sj..] 

A,  ' 

.  l^n]  . 

K.] 

Ae 

where 


(12) 
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(13) 


[•S'zj]  — 


Again,  the  order  of  the  submatrices  [^n]i 

and  [B~]  are  N^xl,  N^xl,  N,xl,  N^xl,  N^xl,  and  NjXl,  respectively.  It 
is  noted  that  the  problem  of  determining  the  matrix  [5i2]  is  related  to  that 
of  [5'2i]  by  reciprocity.  The  explicit  relationship  is  given  in  the  next  section. 

The  scattering  matrix  [522]  is  the  modal  reflection  coefficient  associated 
with  the  interaction  at  the  open  end  of  the  waveguide  which  relates  the 
amplitude  (A“  or  B~)  of  the  n‘^  (TM^  or  TE^)  mode  which  is  reflected 
from  the  open  end  when  a  (TM*  or  TE^)  mode  of  amplitude  (A^  or 
)  is  incident  at  the  open  end.  Thus, 


]  ^  [  («“.]  K,]  1  [  [^,i ' 

.  J  I  Kl  K%]  j  l  kl . 

where 


(14) 


[•S’22]  = 


K1  [O' 
O  O. 


(15) 


Also,  notice  that  the  submatrices  [-Sn'qjj  [-Rnt?]?  are  of 

order  N^xN/,,  Nf,xNe,  N^xN^,  and  N^xN,,  respectively. 
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The  scattering  matrix  [Sp]  is  the  modal  reflection  coefficient  associated 
with  the  termination  at  junction  II  which  relates  the  coefficients  or  am¬ 
plitudes  of  the  waveguide  modes  E~  incident  on  the  termination  to  the 
coefficients  of  the  waveguide  modes  reflected  from  that  termination. 
Thus, 


’  [^:] ' 

'  [c] 

K-]  ] 

.  [rtl 

[ni,l  J 

Ki  J 

(16) 


where 

[rlr,]' 

.[C] 

Again,  the  submatrices  [r**,],  [T",]  are  of  order  N^xN^, 

NhXNe,  NjXNh,  and  N^xN*,  respectively. 

At  any  given  frequency,  the  waveguide  can  support  a  finite  number  of 
propagating  modes  and  an  infinite  number  of  evanescent  modes.  However, 
due  to  reasons  similar  to  those  discussed  in  [9],  only  a  finite  number  of 
evanescent  modes  will  contribute  significantly  to  the  radiation  problem  as 
6  — >  0,  so  that  only  a  finite  number  of  evanescent  modes  needs  to  be 
included  in  practice. 

It  is  observed  from  Figure  4  that  one  may  describe  the  field  scattered 
by  the  cavity  as  consisting  of  two  parts.  The  entire  scattered  field  thus 
consists  of  the  contribution  from  just  the  open  end  by  itself,  together  with 
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the  field  radiated  from  the  open  end  as  a  result  of  all  the  interactions 
between  the  open  end  and  the  termination  that  arise  from  the  incident 
field  which  is  coupled  into  the  waveguide  cavity  region.  Referring  to  Figure 
4,  it  is  observed  that  the  incident  electric  field  E'  at  the  open  end  is  partly 
scattered  by  junction  I  to  produce  a  scattered  field  E*°  given  by  (6)  and 
(8).  The  rest  of  the  incident  field  is  also  coupled  at  I  into  the  waveguide 
region  provided  that  there  is  no  other  scattering  from  the  external  part  of 
the  cavity;  this  field  then  becomes  incident  at  junction  II  from  which  it  is 
subsequently  reflected.  The  field  reflected  from  junction  II  is  incident  back 
at  I  where  it  undergoes  further  scattering  into  the  exterior  and  interior 
regions,  and  so  on,  thus  giving  rise  to  multiple  wave  interactions  between 
junctions  I  and  II.  The  fields  resulting  from  these  multiple  interactions  may 
be  expressed  in  a  convergent  Neumann  series  as  done  by  Pace  and  Mittra 
[11];  however,  an  alternate  procedure  based  on  a  self-consistent  method 
[10]  leads  directly  to  the  same  result.  The  latter  self-consistent  MSM  based 
analysis  will  be  employed  here. 

Let  [F12]  represent  the  total  field  incident  at  I  from  II  after  taking  all 
the  multiple  interactions  into  account.  Likewise,  let  [£^21]  represent  the 
total  field  incident  at  II  from  I  after  taking  all  the  multiple  interactions 
into  consideration.  Then  the  total  scattered  field  in  the  exterior  region 
denoted  by  [F*]  consists  of  a  superposition  of  [E'"]  and  [E”""].  where  [E”*'’] 
is  the  field  scattered  into  the  exterior  region  when  [E^j  is  incident  on  I. 
Therefore, 
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[£•]  =  [£•'’] +  [E^°]  (18) 

where  [£’"*'’]  may  be  expressed  as 


[En  =  l5n][f:n]L.o-  (19) 

The  expressions  for  [  E12  ]  and  [  E21  )  niay  be  written  in  terms  of  the 
matrices  defined  above;  namely, 


|E=.ll„0  =  {iSnKB'l  +  (20) 

and 


=  l^’l|Srl|/’l(£,.|l„„  (21) 

in  which  [Sr]  denotes  the  reflection  coefficient  scattering  matrix  pertaining 
to  junction  II  as  defined  in  (17),  and  [P]  is  a  diagonal  matrix  accounting 
for  the  phase  path  delay  along  the  length  L  between  junctions  I  and  II. 
Eliminating  [£^12]  between  (20)  and  (21)  yields 


(m  -  |S!=l|f||Sr||P|)  IBnII.,.  -  ISnIlP'J.  (22) 

where  [I]  is  an  identity  matrix.  Hence,  it  follows  that 
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(23) 


|£!.)L=o  =  (Ifl  -  ISjiimiSrllPI)-'  |S2,)|Ei. 

From  (21),  becomes 

lEn|l„o  =  l^’llSrl(P|(|/|  -  |5„liJ’||5rl|P))-'  |S„11£')-  (24) 

Incorporating  the  above  result  into  (19)  yields 

|£”"1  =  |S.,l|Pl|5r|lP|(I/|  -  (S:.|lP||Srl|Pl)-'  |S,.]|£’).  (25) 

Finally,  combining  (6)  and  (19)  according  to  (18)  leads  to  the  self-consistent 
expression  for  the  total  scattered  field  [£*]  as 

(£‘l  =  {(5..I  +  |S„|(/>|(5rl|P|((/i  -  |S!!l|fJ|Sr)|P))-'  (26) 

The  scattering  matrices  [^n],  [5i2],  (52i],  and  [522]  can  be  found  in  a 
relatively  simple  manner  via  the  GTD  and  ECM  as  shown  in  the  next 
section. 

Ill  Analysis  for  Determining  the  Elements 
of  the  Scattering  Matrices 

The  elements  of  [5ii]  can  be  found  via  a  direct  application  of  GTD  which 
provides  the  fields  diffracted  from  the  edges  at  the  open  end  of  the  cavity 
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when  a  plane  wave  is  incident  on  these  edges.  Using  the  concept  of  modal 
rays  one  can  also  find  the  elements  of  [Sjj]  directly  using  the  GTD.  On  the 
other  hand,  it  is  necessary  to  employ  GTD  based  ECM  to  find  the  elements 
of  [S21]  and  [522];  however,  [52i]  is  related  to  [5’i2]  via  reciprocity  and  it  can 
therefore  be  found  by  invoking  reciprocity  once  [S12]  is  known. 

For  the  dielectric  loaded  cavity  considered  in  this  report,  the  contri¬ 
bution  to  the  scattered  field  arising  from  the  multiple  wave  interactions 
between  edges  is  significant  outside  the  specular  direction.  Therefore,  the 
field  doubly  diffracted  across  the  open  end  is  also  included  in  the  [5ii], 
[•S'n],  [•522]  calculation  of  the  scattered  field. 

It  has  also  been  found  that  when  the  incident  angle  is  away  from  graz¬ 
ing  the  ground  plane,  a  simple  formulation  for  the  doubly  diffracted  field 
across  the  open  end  can  be  calculated  by  the  equivalent  current  method. 
However,  at  and  near  grazing  incidence,  this  simple  formulation  for  the 
doubly  diffracted  field  is  no  longer  valid.  In  this  case,  a  more  elaborate 
formula  by  Tiberio  and  Kouyoumjian  [12]  (TK)  as  shown  in  Appendix  B 
has  to  be  used  for  calculating  the  doubly  diffracted  field.  The  TK  doubly 
diffracted  field  properly  compensates  the  singularity  of  the  singly  diffracted 
field.  However,  the  TK  formulation  is  much  more  complicated  which  results 
in  a  much  longer  computation  time  and  the  results  are  essentially  the  same 
as  those  obtained  from  the  simple  formulation  mentioned  above  when  the 
incident  angle  is  away  from  grazing  the  ground  plane.  It  is  for  this  reason 
that  the  doubly  diffracted  field  in  this  report  is  given  in  terms  of  the  simple 
GTD  formulation,  and  for  near  grazing  incidence  these  doubly  diffracted 
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fields  can  be  replaced  by  the  TK  formulation  as  in  Appendix  B. 

It  is  noted  that  with  the  TK  formulation,  the  singularity  in  the  singly 
diffracted  field  at  grazing  incidence  can  been  properly  corrected.  But  the 
singularity  in  the  doubly  diffracted  field  still  remains.  This  additional  sin¬ 
gularity  in  the  doubly  diffracted  field  can  ageiin  be  corrected  by  including 
higher  order  (triply,  quadruply,  •••)  diffracted  fields.  But  then,  the  com¬ 
putation  will  become  too  cumbersome  to  be  useful  for  practical  purposes. 
To  avoid  the  need  of  including  the  triply  and  higher  order  diffracted  fields, 
a  cosine  factor  is  introduced  into  the  doubly  diffracted  field  in  the  exterior 
region.  This  cosine  factor  forces  the  doubly  diffracted  fields  to  vanish  along 
the  part  of  the  ground  plane  where  the  doubly  diffracted  fields  alone  do  not 
satisfy  the  boundary  conditions.  Specifically,  this  is  done  as  follows:  Divide 
the  exterior  free  space  into  two  quadrants.  Each  doubly  diffracted  ray  will 
properly  pair  up  with  the  singly  diffracted  ray  and  therefore  no  modification 
is  needed  in  one  of  the  quadrants  that  encloses  the  edge  where  the  doubly 
diffract'*d  ray  emanates.  In  the  other  quadrant,  the  doubly  diffracted  ray 
becomes  singular  at  grazing  incidence  without  proper  compensation  and  so 
a  cosine  factor  cos^,  0  <  <  7r/2,  is  included,  where  ^  =  0°  is  chosen  to 

be  the  boundary  between  the  two  quadrants  and  the  modification  factor  is 
one  at  the  boundary  and  zero  in  the  singular  direction. 

In  addition,  it  is  also  found  that  evanescent  modes  with  a  high  reflection 
coefficient  {5r]  also  contribute  significantly  to  the  scattered  field  outside 
the  specular  region  and  therefore  they  must  be  included  in  the  calculation. 
From  the  cases  treated  here,  it  is  found  that  generally  only  the  first  5  or  6 
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evanescent  modes  are  needed  for  accuracy.  In  order  to  account  for  the  effect 
of  the  evanescent  modes  by  using  the  GTD  and  ECM,  it  has  been  necessary 
to  generzdize  Keller’s  edge  diffraction  coefficient  (see  Equation  (29) 
below)  to  include  complex  angles  of  incidence  and  diffraction. 

III.l  Development  of  [5ii] 

The  matrix  [5n]  describes  the  field  scattered  only  from  the  open  end  of  the 
cavity.  The  elements  of  [5ji]  can  be  found  via  the  GTD  as  indicated  earlier. 
Since  there  is  no  coupling  between  the  z  and  6  polarizations  for  the  2-D 
case,  the  off-diagonal  elements  S*g  and  in  [5ii]  vanish.  The  diagonal 
elements  5'°  and  Sgg  can  be  found  via  GTD  as  follows; 

The  field  scattered  from  just  the  open  end  is  given  by 


+  E*”  (27) 

where  E[\  are  the  fields  diffracted  from  the  top  /  bottom  edge  at  a;  =  ±a/2. 
The  E{°  and  Ej"  are  given  in  terms  of  the  GTD  by 

E;”  ~  E’(Qa.j)  .  E— ^  -K  E‘(Qi.2)  .  E,-—  ;  (28) 

V^l.2  V^2,I 

where  E  is  the  dyadic  form  of  the  2-D  version  of  Keller’s  diffraction  coeffi¬ 
cient,  which  is  given  by 


E  =  D,{<i>',4),n)zz  -  Dh{<i>',4>,n)Ms 
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and  Dd  is  the  dyadic  form  of  the  2-D  version  of  the  double  diffraction  coef¬ 
ficient.  For  incident  angle  away  from  grazing  angle  of  the  ground  plane,  a 
simple  double  diffraction  coefficient  based  on  the  equivalent  current  method 
can  be  applied  and  one  finds  that 


Dd  =  D,{4)\4>i2,n)D['iL,4>ii,4>,n)zz  (29) 

where  (/(jj  is  the  angle  of  the  singly  diffracted  ray  incident  from  edge  1  to 
edge  2  and  represent  the  UTD  [3]  diffraction  coefficients.  Then  from 
(6)  it  is  clear  that: 


5.0  ^  +  e„ir +  e,) 

-  ^„rr  -  e,) 

+ - - - -r - -0.) 

y/C  2  2 

gjfco  5(»intf.-sin«,)-jfc,a  _ 

+ - ^ - D,{Tr  -  0i,-^)D['{k„a,-^,n  +  $,)  (30) 

S»o  ^  + 

+  -  ^.,7r  -  <?,) 

|(»in«,-8inJ,)-jfcoa  _  _ 

+ - r - -,.r  -  «.) 

V  o  2  2 

f(«in«,-8inS,)-j<r,a  ^  ^ 

q. -  - -  e,,-)Dl^(koa,  -,7r  4-  ^,).  (31 ) 

y  a  2  2 

The  first  two  terms  in  Equations  (30)  and  (31)  are  the  fields  associated 
with  rays  singly  diffracted  from  each  of  the  two  edges  at  the  open  end; 
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whereas,  the  third  and  fourth  terms  are  the  fields  associated  with  rays 
that  are  doubly  diffracted  from  those  edges.  The  doubly  diffracted  field 
is  produced  when  the  singly  diffracted  field  from  one  edge  is  incident  at 
the  opposite  edge.  The  D,  h{(f>\4))  in  (30)  and  (31)  represent  Keller’s  GTD 
edge  diffraction  coefficients  [1]  with  n  —  3/2;  the  D,  h  are  valid  only  outside 
the  shadow  boundary  transition  regions.  The  represent  the 

uniform  GTD  (UTD)  edge  diffraction  coefficients  [3]  with  n  =  3/2.  It  is 
noted  that  the  D]\  remeiin  valid  within  the  shadow  boundary  transition 
regions  and  reduce  to  D,  ^  outside  those  regions  [3]. 

III. 2  Development  of  [S12] 

The  radiation  matrix  [Sijj  describes  the  radiation  of  the  waveguide  cavity 
modes  from  the  open  end  into  the  external  free  space.  The  elements  of 
this  matrix  are  found  through  the  use  of  the  GTD;  here,  the  idea  is  to 
decompose  the  waveguide  modes  incident  at  the  open  end  into  modal  rays. 
Each  modal  ray  is  tree  ed  as  a  plane  wave  incident  on  the  appropriate  edge 
from  the  modal  ray  direction  which  is  diffracted  from  the  edge  into  the 
exterior  region  as  shown  in  Figure  5.  The  radiated  electric  field  has  only 
a  i-component  for  the  TM^  modes  and  only  a  ^-component  for  the  TE* 
modes.  Therefore,  Slg^  and  are  zero;  and  the  elements  and 
can  be  found  to  be 

+  e,) 


17 


Z  yJQ.  Z 

(32) 

L  yj  (X  L 

r-e.) 

Z  \J  0,  Z 

(33) 

where 

-  --CD- 

(34) 

It  is  noted  that  the  effect  of  rays  doubly  diffracted  from  the  edges  is  included 
in  the  above  expressions. 

111.3  Reciprocity  Relationship  between  [5i2]  and  [52i] 

As  discussed  in  the  previous  section,  the  elements  of  the  matrix  [5i2]  can  be 
related  to  that  of  [52ij  by  reciprocity  [13j.  For  the  2-D  problem  discussed 
here,  the  relationship  becomes 

IS..1=  (35) 

111.4  Development  of  [522] 

The  elements  of  [522]  arc  the  modal  reflection  coefficients  of  the  mode 
reflected  from  the  open  end  due  to  a  mode  incident  on  the  open  end. 
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These  modal  reflection  coefficients  for  the  open  end  can  be  found  through 
the  help  of  the  ECM  [6]  in  which  the  strengths  of  the  equivalent  currents 
are  found  from  the  GTD.  In  order  to  apply  the  GTD  and  ECM  methods, 
an  incident  waveguide  mode  is  first  decomposed  into  a  pair  of  bouncing 
modal  rays  for  the  parallel  plate  waveguide  (9].  Each  modal  ray  can  be 
treated  as  a  plane  wave  which  is  incident  on  the  edge  in  the  direction  of  the 
modal  ray  angle  as  in  Figure  5.  Only  one  set  of  modal  rays  is  incident  on 
an  edge  to  produce  diffracted  rays.  The  effect  of  the  field  doubly  diffracted 
between  the  edges  is  included  by  treating  the  singly  diffracted  field  from 
one  edge  as  being  the  field  incident  on  the  second  edge  to  produce  the 
doubly  diffracted  field.  The  edge  diffracted  fields  can  be  found  in  terms 
of  the  GTD  edge  diffraction  coefficients  D,,/,;  these  diffracted  fields  can 
be  viewed  as  being  produced  by  a  pair  of  equivalent  line  currents  at  the 
two  edges.  The  equivalent  currents  can  be  chosen  to  radiate  inside  the 
semi-infinite  waveguide  (and  its  geometric  extension  to  form  an  infinite 
waveguide);  the  amplitudes  of  the  modal  fields  excited  by  these  equivalent 
currents  within  the  original  semi-infinite  portion  of  the  waveguide  then 
constitute  the  required  modal  reflection  coefficients  as  shown  in  [14].  For 
propagating  modes,  nir  <  ka  so  the  modal  angle  =  sin“^(^)  as  shown 
in  Figure  5  is  real.  For  evanescent  modes,  tttt  >  ka  so  the  modal  angle  is 
pure  imaginary.  As  a  result  of  the  orthogonality  property  of  the  waveguide 
modes,  the  elements  and  72^^,  can  be  shown  to  vanish,  and  the  elements 
and  are  given  by  [9] 
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Rtt  =  -- 
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(36) 


(37) 


III. 5  Development  of  [Sp] 

The  reflection  coefficient  [5r]  is  the  modal  reflection  coefficient  associated 
with  the  termination.  In  Figure  4(a)  the  cross  section  of  the  waveguide 
cavity  is  uniform  and  the  interior  termination  is  comprised  of  a  planar,  uni¬ 
form  dielectric  layer  which  is  normal  to  the  waveguide  axis  and  backed  by  a 
short  circuit.  Therefore,  from  the  orthogonality  property  of  the  waveguide 
modes,  there  is  no  coupling  between  the  different  waveguide  modes  across 
the  dielectric  interface,  and  at  the  short  circuit.  Accordingly,  the  matrix 
[5r]  becomes  diagonal  for  this  case.  The  elements  of  [Sp]  are  trivially  ob¬ 
tained  by  matching  the  appropriate  boundary  conditions  on  the  tangential 
fields  at  the  dielectric  interface  and  the  short  circuit,  respectively,  and  by 
invoking  the  orthogonality  of  the  modes.  Thus,  it  can  be  shown  that: 
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IV  Numerical  Calculations  and  Discussion 

The  accuracy  of  the  present  hybrid  approach  for  analyzing  the  scattering 
from  the  dielectric  filled  cavity  in  a  ground  plane  is  tested  by  comparing 
the  numerical  results  obtained  using  this  approach  with  those  based  on 
an  independent  numerical  moment  method  (MM)  solution  of  an  integral 
equation  for  the  same  problem.  Several  numerical  results  illustrating  this 
comparison  are  given  in  Figures  6  through  12.  It  is  seen  from  these  figures 
that  the  results  for  the  bistatic  scattering  from  the  cavity  can  be  signifi¬ 
cantly  improved  for  aspects  near  the  ground  plane  by  including  the  doubly 
diffracted  rays  which  result  from  the  interaction  between  the  two  edges  at 
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the  open  end.  The  latter  comment  applies  to  both  types  of  incident  polar¬ 
ization.  For  near  grazing  angles  of  incidence  on  the  cavity,  the  scattered 
field  become?  less  accurate  within  15  degrees  or  so  off  the  ground  plane. 
This  inaccuracy  is  further  enhanced  for  aspects  close  to  the  ground  plane 
when  the  dielectric  constant  is  made  extremely  large  in  the  TM^  case  as 
is  evident  from  Figure  12.  Although  the  use  of  the  TK  double  diffraction 
formulation  generally  yields  better  results  at  near  grazing  incidence  than 
those  obtained  with  the  simple  GTD  double  diffraction  formulation,  it  still 
does  not  completely  remove  the  inaccuracy  problem  at  aspects  close  to  the 
ground  plane.  Possibly,  the  triple  and  other  higher  order  multiple  interac¬ 
tions  could  provide  an  additional  improvement;  however,  the  calculations 
to  include  these  become  more  cumbersome  and  will  not  be  attempted  here. 

Additional  numerical  results  for  the  backscattering  from  the  rectangular 
cavity  in  the  ground  plane  are  also  provided  in  Figures  13-16.  The  latter 
calculations  are  performed  for  different  widths  and  depths  of  the  cavity,  as 
well  as  for  different  dielectric  constants  of  the  material  filling  the  cavity. 
Again,  for  the  aforementioned  reasons,  the  calculations  in  those  figures 
should  remain  accurate  except  to  within  15  degrees  or  so  of  the  ground 
plane. 

In  conclusion,  an  efficient  hybrid  method  combining  the  high  frequency 
techniques  with  the  modal  approach  within  the  framework  of  the  self  con¬ 
sistent  multiple  scattering  formulation  (i.e.  MSM  or  equivalently  GSMT) 
has  been  developed  to  analyze  the  scattering  by  two  dimensional  dielectric 
filled  rectangular  cavities  in  a  ground  plane.  This  approach  can  be  extended 
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to  treat  the  scattering  by  corresponding  three  dimensional  dielectric  filled 
rectangular,  circular  or  other  cavity  shapes  for  which  the  waveguide  cav¬ 
ity  modes  can  be  found  analytically  in  closed  forms.  Furthermore,  it  is 
also  possible  to  extend  this  useful  hybrid  technique  to  include  the  effects 
of  antennas  placed  within  such  cavities  by  combining  it  with  the  moment 
method.  The  latter  topics  are  proposed  as  part  of  a  future  study. 
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Appendix  A 


Waveguide  Modes 


Ct  =  ie,  +  yty 


(43) 


ht  =  x/ij  +  yhy 


A.  Soft  Case  (i.e.,  TM^  Case): 


e,{x,y)  =  A„sin[^(y  + 


a 


.  ,  S  JTITT  ,  r'nTT  ,  a,,  .-a  _ 

^*(a;,y)  =  - ^^nCos[— (y  + 


^rq,nl^n^- 


hy{x,y)  =  ±- — >l„sin[— (y  + 


'eg,n 


e*(a:,y)  =  ey(x,y)  =  /i,(x,y)  =  0. 
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B.  Hard  Case  (i.e.,  TE^  Case): 


Where 


k,(x,y)  =  A„cos[—{y+ 

a  2 

^x{x,y)  =  An  sm\^ {y  + 

p„a  a  2 

f’y(l,y)  =  TZeq.nAnCOs[  —  {y  + 

a  Z 

hAx,y)  =  hy{x,y)  =  e,(x,y)  =  0. 
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Appendix  B 


The  Uniform  GTD  Solution  for  the  Doubly 
Diffracted  Field  between  Two  Edges 

The  Uniform  GTD  (UTD)  solution  for  the  doubly  diffracted  field  between 
two  edges  has  been  found  by  Tiberio  and  Kouyoumjian  [12]  recently.  This 
solution  can  properly  account  for  the  singularity  of  the  singly  diffracted 
field  even  when  the  incident  angle  is  near  the  grazing  angle  of  these  two 
edges,  the  detail  of  this  solution  can  be  found  in  [12].  In  this  appendix,  one 
lists  only  2-D  version  dyadic  form  of  the  solution  wliich  is  given  by 


DaiQuQi)  =  zzDd,  - 


(51) 


where 


Dd,,dh 

=  p..JQi,Q2) 

(52) 

P Ql) 

(53) 

'I 

II 

O' 

II 

&. 

■  PA^P,^J 

(-1)"+’ 

STrjtnim)^  sin*  ^  sin*  ^ 

F(fclap)  -  F(kla,)  6--’''“ 

(54) 

Op  —  Og  y/k^a 

~  — - - -cot  ^  cot  ^ 

4?r;(r?i7j2)  2nj  2^2 
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The  angles  0i, 


I 


a^F{kLap)  —  apF{kLag) 

®p  —  o,  y/ie^a 

(55) 

$ 

=  '^  + p  =  i,4 

(56) 

=  TT  +  (-l)«[<^2  -  9  =  1,4 

(57) 

Op 

=  2cos'[(#p-7r~(-l)'>2n,7r^•p)/2.] 

(58) 

a. 

=  2cos'[($, -7r~(-l)92nj7rJV,)/2.] 

(59) 

Np 

is  the  integer  that  most  nearly  equal8(- 

N, 

is  the  integer  that  most  nearly  equals(- 

1)9  ^9 
2Tl27r 

Vi 


1  for  t  =  1,2 
[  -1  for  i  =  3,4. 

T]i  for  TM,  case,  i.e.  P, 

1  for  TE,  case,  i.e.  P„, 

<f>2  and  <^2  are  defined  in  Figure  17. 
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27 


References 


[1]  J.  B.  Keller,  “Geometrical  Theory  of  Diffraction,”  Journal  of  Optical 
Society  of  America,  Vol.  52,  pp  116-130,  1962. 

[2]  S.  N.  Karp  and  J.  B.  Keller,  “Multiple  Diffraction  by  an  Aperture  in 
a  Hard  Screen,”  Optica  Acta,  Vol.  8,  pp  61-72,  January  1961. 

[3]  R.  G.  Kouyoumjian,  “The  Geometrical  Theory  of  Diffraction  and  its 
Applications,”  in  Numerical  and  Asymptotic  Technique  in  Electromag¬ 
netics^  R.  Mittra,  ed.  New  York,  Springer  Verlag,  1975. 

[4]  P.  H.  Pathak,  “Techniques  for  High  Frequency  Problems,”  Chapter  in 
Handbook  of  Antenna  Theory  and  Design,  eds.  Y.  T.  Lo  and  S.  W. 
Lee,  to  be  published  by  ITT-Howard  W.  Sams  and  Co.,  Inc, 

[5]  R.  Mittra  and  S,  W.  Lee,  Analytical  Techniques  in  the  Theory  of 
Guided  Waves,  The  Macmillan  Company,  New  York,  1971. 

[6]  R.  G.  Kouyoumjian,  P.  H.  Pathak  and  W.  D.  Burnside, “A  Uniform 
GTD  for  the  Diffraction  by  Edges,  Vertices  and  Convex  Surfaces,”  in 
Theoretical  Methods  for  Determining  the  Interaction  of  Electromag¬ 
netic  Wave  with  Structures,  J.  K.  Skwirzynski,  ed.,  Netherlands,  Si- 
jthoff  and  Noordhoff,  1981. 

[7]  W.  D.  Burnside  and  L.  Peters,  Jr.,  “Axial  RCS  of  Finite  Cones  by  the 
Equivalent  Current  Concept  with  Higher  Order  Diffraction,”  Radio 
Science,  Vol.  7,  No.  10,  pp.  943-948,  October,  1972. 


28 


[8]  E.  F.  Knott  and  T.  B.  A.  Senior,  “Comparison  of  Three  High  Fre¬ 
quency  Diffraction  Techniques,”  Proceeding  of  IEEE,  Vol.  62,  pp.  1468- 
1474,  1974. 

[9]  A.  Altmtas,  “Electromagnetic  Scattering  from  a  Class  of  Open-Ended 
Waveguide  Discontinuities,”  Ph.  D.  Dissertation,  The  Ohio  State  Uni¬ 
versity,  Department  of  Electrical  Engineering,  Columbus,  Ohio,  Win¬ 
ter  1986. 

[10]  C.  C.  Huang,  P.  H.  Pathak,  C.  Y.  Lai  and  D.  L.  Moffatt,  “Analysis 
of  Electromagnetic  Backscatter  from  an  Inlet  Cavity  Configuration,” 
Final  Report  712661-4,  October  1982,  The  Ohio  State  University  Elec- 
troScience  Laboratory,  Department  of  Electrical  Engineering;  prepared 
under  Contract  No.  F19628-80-C-0056  for  the  Department  of  the  Air 
Force,  Electronic  Systems  Command,  Hanscom  Air  Force  Base. 

[11]  J.  Pace  and  R.  Mittra,  “Generalized  Scattering  Matrix  Analysis  of 
Waveguide  Discontinuity  Problems,”  Quasi-Optics,  Vol.  XIV,  pp.  177- 
197,  Polytechnic  Institute  of  Brooklyn  Press,  New  York,  1964. 

[12]  R.  Tiberio,  G.  Manara,  G.  Pelosi  and  R.  G.  Kouyoumjian,  “High- 
Frequency  Diffraction  by  a  Double  Wedge,”  to  be  published  in  IEEE 
APS. 

[13]  P.  H.  Pathak,  C.  W.  Chuang  and  M.  C.  Liang, “Inlet  Modeling  Stud¬ 
ies,”  Technical  Report  717674-1,  October  1986,  The  Ohio  State  Univer¬ 
sity  Electroscience  Laboratory,  Department  of  Electrical  Engineering; 


29 


prepared  under  Contract  No.  N60530-85-C-0249  for  the  Department 
of  the  Navy,  Navy  Weapons  Center,  China  Lake,  CA  93555. 

[14]  P.  H.  Pathak  and  A.  Altintas,  “An  efficient  High  Frequency  Analy¬ 
sis  of  Modal  Reflection  and  Transmission  Coefficients  for  a  Class  of 
Waveguide  Discontinuities,”  to  appear  in  Radio  Science. 


30 


CO  CO 
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Figure  3:  Isolated  interactions  at  junctions  I,  II  and  III. 
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Figure  4:  The  Interactions  of  the  Multiple  Scattering  Matrix. 
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-90.  -60.  -30.  0.  30.  60,  90. 


(a)  TE,  Case 

—  MM, - AMSM  (with  Simple  Double  Diff.), 

-  •  -  AMSM  (with  Single  Diff.  only)  and 
-  -  AMSM(  with  Modif.  TK  Double  Diff.). 

Figure  6:  Comparison  Between  Moment  Method  and  AMSM  6- Approach 
for  the  Bistatic  Scattering  Case  with  a  =  2.8  A,  L  =  0.25  A,  e,.  =  2.49,  /x, 
=  1  and  0i  =  0°. 
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(b)  TM,  Case 


Figure  6:  (Continued) 


37 


-60.  -30.  0. 


30.  60. 


-60.  -30.  0.  30.  60.  ? 
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(a)  TEj  Case 

—  MM,  -  -  -  AMSM  (with  Simple  Double  Diff.), 

-  •  -  AMSM  (with  Single  Diff.  only)  and 
—  AMSM  (with  Modif.  TK  double  Diff.). 


Figure  7:  Comparison  Between  Moment  Method  and  AMSM  ^-Approach 
for  the  Bistatic  Scattering  Case  with  a  =  2.8  A,  L  =  0.25  A,  e,  =  2.49,  /Xr 
=  1  and  9,  =  -30°. 
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(b)  TMj  Case 
Figure  7:  (Continued) 
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Figure  8:  Comparison  Between  Moment  Method  and  AMSM  6-Approach 
for  the  Bistatic  Scattering  Case  with  a  =  2.8  A,  L  =  0.25  A,  =  2.49, 

=  1  and  6i  =  -60°.  40 
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(b)  TM*  Case 
Figure  8:  (Continued) 
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'-90.  -60.  -30.  O.  30.  60.  90. 

ANGLE  IN  DEGREES 

(a)  TEj  Case 

—  MM, - AMSM  (with  Modif.  Simple  Double  Diff.)  and 

-  •  -  AMSM  (with  Modif.  TK  Double  Diff.  ). 

Figure  9:  Comparison  Between  Moment  Method  and  AMSM  6-Approach 
for  the  Bistatic  Scattering  Case  with  a  =  2.8  A,  L  =  0.25  A,  e,  =  2.49,  p, 
=  1  and  0t  =  -75°.  42 
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Figure  9:  (Continued) 
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-90.  -60.  -30.  0.  30.  60.  90. 


(a)  TE,  Case 

—  MM, - AMSM  (with  Modif.  Simple  Double  Diff.) 

and  -  •  -  AMSM  (with  Modif.  TK  Double  Diff.). 


Figure  10:  Comparison  Between  Moment  Method  and  AMSM  6-Approach 
for  the  Bistatic  Scattering  Case  with  a  =  2.8  A,  L  =  0.25  A,  e,  =  2.4.9,  fir 
=  1  and  0i  =  -80°.  44 
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Figure  10:  (Continued) 
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ANGLE  IN  DEGREES 

(a)  TEj  Case 

—  MM, - AMSM  (with  Modif.  Simple  Double  Diff.) 

-  ■  -  AMSM  (with  Modif.  TK  Double  Diff.) 

Figure  11:  Comparison  Between  Moment  Method  and  AMSM  for  a  very 
shallow  rectangular  cavity  with  L  =  0.15  A  a  =  2.8  A.  The  Calculation  is 
for  the  Bistatic  Scattering  Cfise  with  =  1,  c,  =  2.49  and  6i  =  -60°. 


Figure  11:  (Continued) 
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-  •  -  AMSM  (with  Modif.  TK  Double  Diff.) 

Figure  12:  Comparison  Between  Moment  Method  and  AMSM  for  a  high 
dielectric  loading.  The  Calculation  is  for  the  Bistatic  Scattering  Case  with 
a  =  2.8  A,  L  =  0.25  A,  =  1  and  =  0°. 
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(b)  TMj  Case  with  e,  =  60 
Figure  12:  (Continued) 
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Figure  12:  (Continued) 
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(a)  TEi  Case 

Figure  13:  The  AMSM  ^-Approach  Calculation  for  the  Backscatter  with  a 
=  2.8  A,  L  =  0.75  A,  e,  =  2.49,  and  Hr  =  1. 
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Figure  13:  (Continued) 
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-75.-60.-45.-30.-15.  0.  15.  30.  45.  60.  75. 


Figure  14:  The  AMSM  ^-Approach  Calculation  for  the  Backscatter  with  a 
=  4.2  A,  L  =  0.25  A,  e,  =  2.49,  and  p,  =  1. 


53 


75.-60.-45.-30.-15.  0.  15.  30.  45.  60. 

ANGLE  IN  DEGREES 

(b)  TM,  Case 
Figure  14:  (Continued) 
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Figure  15:  The  AMSM  6-Approach  Calculation  for  the  Backscatter  with  a 
=  2.8  A,  L  =  0.25  A,  c,  =  9.,  and  /x,  =  1. 
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Figure  15:  (Continued) 
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Figure  16:  The  AMSM  ^-Approach  Calculation  for  the  Backscatter  with  a 
=  4.2  A,  L  =  0.75  A,  =  9.,  and  fir  =  1. 
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Figure  17:  The  Geometry  for  the  Doubly  Diffracted  Field  Between  Two 
Edges 


59 


